Feynman's Propagator Applied to Network Models of Localization 
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Network models of dirty electronic systems are mapped onto an interacting field theory of lower 
dimensionality by interpreting one space dimension as time. This is accomplished via Feynman's 
interpretation of anti-particles as particles moving backwards in time. The method developed maps 
calculation of the moments of the Landauer conductance onto calculation of correlation functions of 
an interacting field theory of bosons and fermions. The resulting field theories are supersymmetric 
and closely related to the supersymmetric spin-chain representations of network models recently 
discussed by various authors. As an application of the method, the two-edge Chalker-Coddington 
model is shown to be Anderson localized, and a derealization transition in a related two-edge 
network model (recently discussed by Balents and Fisher) is studied by calculation of the average 
Landauer conductance. 
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I. INTRODUCTION 

Dirty electronic systems exhibit a variety of interest- 
ing phases and transitions in their transport properties 
at low temperature. For example consider two dimen- 
sional electrons moving in a random potential and a 
strong perpendicular field, as in a quantum Hall experi- 
ment. Generically the wavefunctions are localized, with 
tails that decay over a length scale called the localization 
length, which leads to insulating behaviour at low tem- 
perature. However if the parameters (e.g., the magnetic 
field or electron density) are tuned to special isolated val- 
ues the localization length diverges. Such derealization 
transitions were postulated shortly after the discovery 
of the quantum Hall effect and invoked in order to ex- 
plain it |lj . Since then impressive strides have been made 
in experimental characterization of these transitions ; 
but a satisfactory theoretical understanding has not yet 
emerged ||. 

An important theoretical advance was made by 
Chalker and Coddington who introduced a simplified net- 
work model of the quantum Hall transition 0] . Numeri- 
cal studies show this model produces the same universal 
behaviour at the derealization transition as more literal 
(and more complicated) models of the quantum Hall sys- 
tem p|. Because the network model is relatively simple 
and is based on a clear physical picture of the transi- 
tion, it seems a promising starting point for a controlled 
approximate analysis of the transition. 

More recently it has been observed that the random 
bond Ising model is closely related to a variation on the 
network model |^| . Progress in analysis of network mod- 
els is therefore desirable from this point of view also. 

Conductance is a sensitive probe of derealization. The 
purpose of this paper is to introduce a technique suit- 
able for calculation of the conductance of network mod- 
els. Following Landauer we imagine electrons are injected 
from the source into the sample where they undergo mul- 



tiple scattering |6| . Eventually an electron may either get 
scattered forward into the drain or it may get scattered 
backward into the source (see fig 1). We wish to calculate 
the probability of forward scattering. 



It is often fruitful in statistical mechanics to map a 
problem onto a quantum field theory of lower dimension- 
ality by reinterpreting one space dimension as time. In 
applying that strategy here it becomes neccessary to take 
into account the fact that electrons will then appear to 
move both forwards and backwards in "time" . 
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> 

x-axis 



Figure 1. Landauer's picture: The electron undergoes 
multiple scattering in the sample and is eventually scat- 
tered forward into the drain. In the field theory rep- 
resentation this corresponds to a process involving pair 
creation and annhilation. 
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A similar situation is encountered in quantum elec- 
trodynamics: in a famous paper, Theory of Positrons, 
Feynman showed that it is possible to regard positrons 
as electrons moving backwards in time 0. Following 
Feynman it is possible here to interpret the x-axis in fig 
1 as time provided we regard an electron moving to the 
left as a sort of anti-particle of an electron moving to 
the right. For example the process in which a right mov- 
ing electron is scattered to the left could be regarded as a 
process in which the particle meets it's anti-particle lead- 
ing to their mutual annhilation. From this perspective, 
the simple process depicted in fig 1 where an electron 
zig-zags through the sample could be regarded as a pro- 
cess in which a particle-antiparticle pair is created (at 
the earlier kink in the trajectory) and annhilated (at the 
later kink). 



Based on this interpretation, it is possible to map a 
two-dimensional network model onto a one-dimensional 
field theory of particles and anti-particles. The statistics 
of the particles may be taken to be either Bose or Fermi. 
For technical reasons it is most convenient to introduce 
both species. Thus calculation of the conductance is 
mapped onto calculation of the correlation functions of 
an interacting field theory of bosons and fermions. 



The second model was recently discussed in context of 
quantum transport by Balents and M.P.A. Fisher p3| . It 
is related to a model of glass first studied by Dyson |TJ] 
and to a special version of the random bond Ising model 
introduced by McCoy and Wu jl5) . An enlightening dis- 
cussion of these connections, with references, is given in 
the paper of Balents and M.P.A. Fisher. This model is 
known to have a critical point at a special value of its 
parameters which is of considerable interest as a sim- 
ple example of a random quantum critical point. Much 
is known about this model and its relations, particu- 
larly through the application of real space renormaliza- 
tion group methods by D.S. Fisher jig . In their recent 
paper Balents and M.P.A. Fisher have applied supersym- 
metry methods to calculate the exact two-parameter scal- 
ing function of the Green's function of this model. Here 
we shall study it's derealization transition by calculat- 
ing the conductance. A summary of our results for this 
model is given in section IVC. 

II. FIELD THEORY FORMULATION 

In this section a simple two-edge network model is in- 
troduced and it's Landauer conductance defined. The 
main result is an exact formal expression for the Lan- 
dauer conductance, eq (42-44), within a field theory for- 
mulation of the problem. 



A. Landauer Conductance 



For clarity we illustrate the method on two-edge net- 
work models which are essentially one-dimensional. They 
map onto quantum mechanics problems of interacting 
bosons and fermions — a zero-dimensional field theory. 
Generalization of the mapping to a two-dimensional 
network model is straightforward. Indeed the method 
should be much more broadly applicable to dirty elec- 
tronic systems. Since an important obstacle to non- 
perturbative analysis of random systems has been the 
lack of suitable representations of the problem, it is hoped 
this method may prove useful. The present method is 
closely related to supersymmetric spin-chain representa- 
tions of network models that have recently been discussed 
by several authors Jq-pJ. 



As an application of the method we analyse two dif- 
ferent network models. The first is essentially a one- 
dimensional Chalker-Coddington model which exhibits 
Anderson localization: for a large enough sample the zero 
temperature conductance is found to decay exponentially 
with sample size. 



The models considered in this paper consist of two 
counter-propagating "edge states" coupled by tunnelling. 
Along each one-dimensional edge electrons can propa- 
gate in only one direction (see fig 2). The electron wave 
function has two components: ip+(x), the amplitude to 
be on the right moving edge at position x, and ip-(x), 
the amplitude to be on the left moving edge at position 
x. The time independent Schrodinger equation governing 
this model is 



m(x) 
m*{x) i£ 



= E 



(1) 



The tunnelling amplitude m(x) is some given function. 
Eventually we will take tunnelling to be a random pro- 
cess and will be interested in performing averages over an 
ensemble of different realizations of m{x) with statistics 
to be described below. We are interested in solutions at 
some fixed energy E (the fcrmi energy). 

For later reference, it is useful to rewrite eq (I) in the 
form 



' dx V V J 



E 
m*{x) 



-m{x) 
-E 



(2) 
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For a given realization of m(x) the connection between 
(7, S) and (a, /?) can be found by integrating eq (2) across 
the sample 



T 



(4) 



T is a 2x2 matrix called the transfer matrix of the sam- 
ple. The transfer matrix is the focus of numerical studies 
of network models |l . 

Here we are interested in solving eq (2) subject to scat- 
tering boundary conditions in the probes so that 



Figure 2 (a) A two edge network model, (b) Schematic 
picture of the random function m(x). 

An obvious generalization of this model is to consider 
2A^ counter-propagating edge states. With suitable sta- 
tistical assumptions about the tunnelling these general- 
izations lead to the Chalker-Coddington model of the 
quantum Hall transition [Q or the random bond Ising 
model ||. The model studied here can be considered a 
special anisotropic case in which tunnelling between al- 
ternate pairs of edges has been turned off. 

Following Landauer's method ||, the model is sepa- 
rated into the sample (region between Xi and xf) and 
the probes (regions to the left of Xi and to the right of 
Xf); see fig 2. The left probe is called the source; the 
right probe, the drain. In contrast to the sample the 
probes are assumed to be disorder free. m(x) is there- 
fore random only for Xi < x < Xf. It is constant in the 
probes. For simplicity we shall assume that the edges 
are uncoupled in the probes (that is, m(x) = 0) 1 . This 
is shown schematically in fig 2 (b). 

The wave eq (2) is now trivially soluble in the probes. 
We find 

/ ip+ \ ( a ex P iEx \ 
I ?/>_ J \ f3 exp— iEx J 



a 1, /3 -► r; 7 -► t, S -► 0. (5) 

According to Landauer the zero-temperature conduc- 
tance of the sample (at Fermi energy E) is given by 




(6) 



Returning to eq (2) it is useful to regard x — > time and 
to interpret it as the time- dependent Schrodinger equa- 
tion of a two-level system. An awkward feature of such a 
reinterpretation is that the 2x2 "Hamiltonian" matrix in 
eq (2) is not Hermitean 2 . Sadly this feature will persist 
through much of our analysis and appears to be more 
generally pervasive ||,[l^,[l3| • 

We can now imagine calculating the retarded prop- 
agator for eq (2). Given the wave function (values of 
tp+ and at a particular a;— slice, the retarded prop- 
agator gives the wave function for "future" values of x. 
For example, if the wave function is known at x%, the 
retarded propagator gives the wave function inside the 
sample (and beyond, in the drain). This is in very much 
the same spirit as the transfer matrix (indeed the re- 
tarded propagator is the transfer matrix for a special 
value of it's arguments). It is not helpful for calculating 
the scattering amplitudes. 



1 It is not essential for our approach to assume m(x) = in 
the probes. However, it simplifies the analysis. It would be 
expected that the transport properties of the sample should 
not be modified by the nature of the probes. This can be 
checked explicitly under some circumstances. 

2 Perhaps this is a good place to emphasize that in this paper 
although the "Hamiltonian" is frequently non-Hermitean, it is 
always a well defined operator in a well defined Hilbert space. 
There are no difficulties of interpretation as our "Hamilto- 
nian" is merely used as a calculational device. The true 
Schrodinger equation of the model, eq (1), is governed by 
a Hermitean Hamiltonian. 
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While solving eq (2) subject to scattering boundary 
conditions we do not know the complete wave function 
for any x— slice. Instead we know the positive frequency 
component V>+ = 1 at the earlier slice X4 and we know the 
negative frequency component ip- — at a later slice x / . 
We would like to reconstruct the wave function between 
the slices using this information; or at least we would like 
to know ip+(x — *■ Xf) = t. 

In Theory of Positrons Feynman showed precisely this 
would be achieved by a modification of the retarded prop- 
agator, now known as the Feynman propagator Feyn- 
man derived an integral equation obeyed by his propa- 
gator and a perturbation series for it; this was the main 
focus of his paper. Peripherally, in two Appendixes, he 
showed the same results could be deduced from a sec- 
ond quantized theory of particles and anti-particles. This 
development will now be used to provide a field theory 
representation of the two edge network model (eq 2) 3 . 



B. Deduction From Second Quantization 



The superscript R is superfluous for the moment; its 
function will become apparent later. In terms of the c R 
fermions the Hamiltonian is 



c^c R ) 



m(x)c 



+ c - 



RR 



+ m*(x)c_c 



(9) 



The scattering amplitude t can be computed straight- 
forwardly using this Hamiltonian; see eq (19) below. In 
the remainder of this subsubsection that result will be 
derived. 

Note that Hp is non-Hermitean — this is traceable to 
the non-Hermiticity of the "Hamiltonian" of eq (2). The 
reader troubled by this non-Hermiticity should again 
read footnote 2. Note also that Hp does not conserve 
the total number of c— fermions. Instead we may regard 
the c_fermion as the "anti-particle" of the c+. What is 
conserved then is the total "charge" — the difference in 
the number of particles and anti-particles, given by 



>f = c 



+ 



R] R. 
c_c_. 



(10) 



The purpose of this section is to derive a second quan- 
tized Hamiltonian from which one can calculate the Feyn- 
man propagator and scattering amplitude t of eq (2). 



1. Fermion Representation 



Return to the interpretation of eq (2) as a time de- 
Introduce a+ which cre- 



pendent two level system, 
ates the right moving state 



which creates 



the left moving state 



; and a+ and a_, the cor- 



responding destruction operators. The a's obey fermi 
(anti-)commutation rules. Also introduce the "time"- 
dependent Hamiltonian 



H{x) = (a^al) 



E 
m* (x) 



-m(x) 
-E 



(7) 



This Hamiltonian is not useful for our purpose. It can be 
shown to generate the retarded propagator. 

Instead we must now introduce c R fermions, related to 
the a fermions via the particle-hole transformation 



c_ = a_ 



(8) 



The S-matrix corresponding to H R is the solution to 

(11) 



i-—Sp(x,Xi) = Hp(x)Sp(x,Xi) 



dx 

subject to the initial condition S R (xi,Xi) = 1. 

Turning to the derivation of eq (19), we first obtain 
a useful formula, following Feynman. Let (e+(x), e_ (x)) 
be a solution to the wave eq (2). Construct the opera- 
tor F = e + (xi)c R ^ ! + e-(x i )c R _. F is transformed by the 
S— matrix as follows: 

S§(x,x i )Psf!;- 1 (x,x i ) = e+(x)cf + e_(x)c R . (12) 

This is Feynman's formula. To prove it, regard eq (12) 
as an ansatz; in other words assume that SFS^ 1 is of 
this form with e+(x) and e-(x) some suitable functions. 
Upon differentiation with respect to i, it will be seen 
that the ansatz is consistent provided (e+,e_) obey the 
wave eq (2). 

Carrying out this plan, we note 



i ^ Sp (x, Xi 



-S R -\x, Xl )H R {x) 



(13) 



obtained by differentiating Sp Sp — 1 and using the 
evolution eq (11) for the S— matrix. Eq (11), (12) and 



(13) together yield 



3 The next two subsections follow Feynman very closely. For 
this reason we have given them the same titles as the corre- 
sponding appendixes in Feynman's paper 
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^(sgFS*- 1 ) =[H Fl S§FS*- 1 ] 



e + (x)[H F ,cf}+e-(x)[H F ,c R }. 

(14) 



The relevant commutators are 



[H F ,c*}] = Ec** +m*d± 
[H F ,c R ] 



Ec R - mcf 



(15) 



Hence the derivative of the left side of the ansatz eq (12) 



is 



— I 



2. Boson Representation 

Tracing through the derivation of eq (19) it is seen 
that the Fermi statistics of the c R particles plays no cru- 
cial role; only the commutators of eq (15) are essential. 
Hence we may replace the c R fermions with b R bosons 
governed by the Hamiltonian 

H R (x) = E{bfb R + b^b R ) + m(x)b^b^ + m*(x)b R b R . 

(20) 

Apart from the replacement c R — > b R , this differs from eq 
(9) in the sign of the pair creation term (the term propor- 
tional to 6**6**). This sign change ensures commutation 



dx 



( r - R—i\ m uiunai id u + u_ ). _liiis sign ciiangc ensures com 

[S F FS F ) = c + [Ee+{x) - m{x)e-(x)} relations of the desired form similar to eq (15): 



+c R [m*(x)e+(x) -Ee-(x)]. (16) 



On comparing with the corresponding a;— derivative of 
the right hand side of the ansatz it is seen that indeed 
(e + ,e_) must obey the wave eq (2); this completes the 
proof of Feynman's formula eq (12). 

Next let (e + ,e_) be the special solution that obeys 
scattering boundary conditions: 



[H§,b*p] =Eb Rt + m*(x)b R , 
[H R ,b R] ] = -Eb R - m(x)b R] . 



(21) 



e+(xi) = 1; e-(xi) = r; 
e+(x f ) =t; e-(x f ) = 0. 

This leads to the scattering formula 

S R (x f , Xi ) (cf + rc R ) S R - 1 (x f , x t ) - tc Rf 

=> S R (x f .Xi)c^ = tc Rf S R (x f ,x i )-rS R (x f ,x i )c R , 



(17) 



Note that due to the sign change in the pair creation 
term the boson Hamiltonian, eq (20), is actually Her- 
mitean. Again the 6_ boson may be regarded as the 
anti-particle of the b + . The Hamiltonian conserves the 
difference in the number of particles and anti-particles, 
that is, the charge 

b R Jb R - b m b R . 



Qi 



(22) 



The entire analysis of the preceding subsubsection can 
now be carried over essentially unchanged. The S-matrix 
corresponding to Hg is the solution to 



g 

^~Q^.^B ( x i x i) — H B (x)S B (x,Xi) 



(23) 



(18) 



a useful special case of Feynman's formula, eq (12). 

Finally analyse the matrix element < 0\c R Sp(xf, Xi)c R ' |0 
here |0 > is the vacuum for the c— fermions. Use of the 
scattering formula, eq (18), reveals 



t = 



<0\c R S R {x f ,Xi)cf\0> 
<0\S R (x f , Xl )\0> 



(19) 



Eq (19) shows that the scattering amplitude t can be 
computed by studying the evolution of the single parti- 
cle state c R '\0 > under the Hamiltonian of eq (9). Note 
that the vacuum amplitude in the denominator of eq (19) 
is not trivial. In Feynman's words, "It differs from unity 
because, for example, a pair could be created which even- 
tually annhilates itself" . 

Eq (19) is the main result of this subsubsection. By 
itself, eq (19) is not especially helpful. Our eventual pur- 
pose is to compute averages over the random tunnelling 
process m(x). Since m(x) appears in both numerator 
and denominator of eq (19), this form is not particularly 
well adapted for averaging. 



subject to the initial condition S R (xi,Xi) = 1. Feyn- 
man's formula, eq (12), and the scattering formula, eq 
(18), apply upon making the replacement c R b R . Us- 
-jjag the scattering formula we can deduce 

f = <0\b R S R (x fl x t )bf\0> 
<0|S|(a: / ,a: i )|0> 
— the bosonic analogue of eq (19). 

C. Analysis of the Vacuum Problem 

Remarkably, the vacuum amplitudes for the bosons 
and fermions cancel; that is, 

< 0\S$(x f ,Xi)\0 >< 0\S§(x f ,Xi)\0 >= 1. (25) 

We have previously derived a fermionic and bosonic ex- 
pression for the scattering amplitude, neither well suited 
for performing an average over the random tunnelling 
process, m(x). Eq (25) will enable us to weld these ex- 
pressions into a form suitable for averaging. We return to 
this point in section E below. Here we focus on proving 
eq (25). 
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To calculate the vacuum amplitude, following Feyn- 
man, let us analyse a series of problems that interpolate 
smoothly between a soluble limit and the problem we 
want to solve. Introduce a truncated problem for which 
m(x) is left unchanged for xo < x < Xf and is set equal 
to zero for Xi < x < x . By varying xq we obtain the 
desired series of problems. Evidently xo = Xi is the case 
we wish to solve. On the other hand, if xq = Xf, there is 
no tunnelling and the problem is trivially soluble. 

Denote the scattering coefficients for the truncated 
problem r(xo) and i(xo). For the soluble case, r(xf) = 
and t(xf) = 1. 

The fermion S-matrix for the truncated problem will 
be written as S^f(x,xo)- It obeys 



d_ 

dx* 



- i^S* F (x,x ) = H*(x)S* F (x,x ) 



(26) 



with S^ F (xo, xo) = 1. The truncated boson S-matrix is 
similarly defined. 

The vacuum amplitude for the truncated fermion prob- 
lem will be denoted CJ?(xq). Thus 



C*(x o )=<0\S% r (x f ,x o )\0> 



(27) 



Analogous expressions can be written for the bosonic 
case. When xo — Xf the vacuum amplitude is unity 
for both bosons and fermions. 

The strategy now is to study the variation of C(x ) 
with xq by writing a differential equation for it. A for- 
mal solution of the differential equation can be obtained 
which will suffice to prove the cancellation of the vacuum 
amplitude. 

For this purpose it is useful to write the solution to eq 
(26) as a formal series 



S* F (x f ,x ) = l + i 



f ' 

i I dx\Hp(xi) 

J x 

+^ 2 f ' dxt f 1 dx 2 H%{x 1 )H%(x 2 ) 

J Xn J Xn 



+ ... (28) 
By differentiation of this series it follows 

- i-^- Q Sg F (x f ,x ) = -S^(x f ,x )HP(x ). (29) 

Taking the vacuum expectation of eq (29) obtain 

-^C™ = < 0\S« F (x f ,x )HK(x o )\0 >; 

Hp(x ) = -m(x )c^ t c Kt + others (30) 

In eq (30) we have explicitly written only that term of 
Hp (eq 9) which does not annhilate the vacuum. 
According to the scattering formula, eq (18), 



S* F (x f ,x )c? = t(x )c?S* F (x f ,x ) 

-r(x )S* F (x f ,x )cK. (31) 

Applying this to the matrix clement in eq (30) yields the 
differential equation obeyed by Cj?(xo): 

- i^-C§(x ) = -m(x )r(x Q )C^(x Q ). (32) 
clxo 

We seek a solution subject to Cp(xj) = 1. 
Define L(x ) as the solution of 



d 

i-^-L(xo) = -m(x )r(x ) 



(33) 



with the initial condition L(xf) — 0. Given m(x) one can 
imagine solving the truncated problem for various values 
of xo and, in principle, computing L(xo). The solution 
of eq (32) is then 



Cp{x ) = expL(x ). 



(34) 



Note that this solution obeys the differential equation 
(32) and the initial condition Cp(xf) = 1. 

The same analysis leads to the boson vacuum ampli- 
tude 



c b( x o) = cxp-L(x ). 



(35) 



The sign change is traceable to the sign difference in the 
pair creation terms of the Fermi and Bose Hamiltonians. 
Hence Cp(xo)C§ (xo) = 1 for all xo and in particular for 
x = Xi, which establishes the desired result, eq (25). 

It is interesting that Feynman appears to have devoted 
some effort to interpreting Dirac particles as bosons even 
though it violated the spin-statistics theorem. In partic- 
ular, the result of eq (25), which is of no obvious value 
in quantum electrodynamics, is derived in the Theory of 
Positrons. 



D. Conjugate Amplitude 

The Landauer conductance is given by \t\ 2 . Thus ex- 
pressions for t* analogous to eq (19) and (24) will be 
needed. To obtain them consider the equation conjugate 
to eq (2): 



.d_ 

dx 



-E m*(x) 
-m(x) E 



0+ 



(36) 



Eq (36) is constructed to have the property that if 
(e + ,e_) is a solution to eq (2), then the complex con- 
jugate (e* + ,e*_) is a solution to eq (36). Hence if eq (36) 
is solved subject to scattering boundary conditions, the 
scattering amplitude will be t*. 

On the other hand, direct comparison reveals that 
eq (2) is transformed into eq (36) by the replacements 
E — > —E,m — > — m*,m* — > — m. Hence introduce c A 
fermions governed by the Hamiltonian 
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H A (x) = -E(cfc A + c^c A ) + m*(x)c?c? - m(x)c A c A . 

(37) 

This Hamiltonian is obtained from eq (9) by making the 
replacements indicated above. By the reasoning that lead 
from eq (9) to (19), Hp generates the conjugate ampli- 
tude via 



<0\c A S A (x f ,x,)cf\0> 
< 0\S$(x f ,Xi)\0 > 



t* 



(38) 



where Sp(xf,Xi) is the S-matrix corresponding to H 
Similarly, one can write 



A 
F ■ 



t* 



<0\b A S A (x f ,x t )bf\0> 



(39) 



<0\S^{x f ,x i )\0> 
where Sg(xf,Xi) is the S-matrix of the boson Hamilto 



man 



H A = -E(b Ai b A + b A h A ) - m*{x)b A] b M - m(x)b A b A . 

(40) 

Finally note that the vacuum amplitudes for the A 
bosons and fcrmions cancel, as they do for their R coun- 
terparts: 

< 0\S A {x f ,Xi)\0 >< 0\S£(x f ,Xi)\0 >= 1. (41) 



E. Supersymmetry 

The results of the previous subsections can now be as- 
sembled into an expression for the Landauer conductance 
suitable for averaging over the random tunnelling process 
m(x). Simultaneously introduce c R and c A fermions and 
b R and b A bosons governed by the total Hamiltonian 



Hsusy(x) = H r + H a + H r + H a 
= EM + m(x)A + m*(x)B 



(42) 



Here M = (cfcf + c m c R -R^ A) + (b Rf b R + b m b R - 
R -> A) and A = -cfc^ - c A c A + bfb** - b A b A and 
B = c R c R + cf c At + b R b R - 6f 6 At . 

As usual the corresponding S-matrix obeys 



- ig^SsuSY(x,Xi) = H SU sY{x)S SUS Y(x,Xi) (43) 

and the initial condition SsusY(xi, Xi) = 1. 

The Landauer conductance can be calculated using 



\t\ 2 =< 0\<4c%SsusY(xt,Xi)<$4*\ > 



(44) 



To see this note that since the different particle species 
don't interact, Ssusy = S R SpS§Sg\ hence the matrix 
element in eq (44) decouples into 



< 0\C R Spcf\0 >< 0\c A S^c Af \0 >< 0|S||0 >< 0|5^|0 > 

Using eq (19), (38), (25) and (41), this product is easily 
seen to be \t\ 2 . 

Eq (44) is our principal tool to analyse the network 
model. It shows that the conductance can be calculated 
by following the evolution of a two fermion state under 
the Hamiltonian Hsusy ■ The principal advantage of this 
expression is due to the absence of a denominator which 
makes it well suited for averaging over the random tun- 
nelling process m[x). The averaging will be discussed 
further in sections III and IV assuming different distri- 
butions for m(x). The reader has perhaps noticed that 
one could develop various other expressions for the Lan- 
dauer conductance that involve evolving, instead of a two 
fermion state, a two boson state or a one boson, one 
fermion state; but eq (44) is the form we shall use here. 

Finally let us briefly examine some symmetries of 
Hsusy- Evidently, it conserves the R fermion charge; 
that is, 



[Qp, Hsusy] 







(45) 



where Q R = c R — c 1 ^ c R . Similarly it also conserves 
the R boson charge as well as the corresponding quantity 
for the A boson and fermion. Hsusy also possesses a su- 
persymmetry analogous to that displayed by the models 
studied in refs [p|-p^[. Specifically the Hamiltonian, eq 
(42), commutes with the supercharges Q r b ,Qp B ,Qbf 
and Qgp which are given by 



Qfb 
Qbf 



cz'b_, 



fb R - cV*; 



(46) 



Qp B and Q BF are obtained by replacing R — > A in eq 
(46). 

It will be seen below that after disorder averaging we 
must analyse an effective interacting Hamiltonian instead 
of the non-interacting (but random) Hamiltonian of eq 
(42). The effective interacting Hamiltonian will also pos- 
sess the symmetries discussed above. 



III. TWO EDGE CHALKER-CODDINGTON 
MODEL 



A. Disorder Average 



First let us assume that the real and imaginary parts 
of m(x) are independent gaussian white noise processes 
with zero mean. Thus [m(x)] cns = 0. Here [. . .] cns de- 
notes an average over the ensemble of m{x) realizations. 
The variance is given by 
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[m*(x), m(x')]ons = D6(x — x') 
[m(x),m(x')] ens = 0. 



(47) 



Essentially the same statistics result if it is assumed 
that m(x) is a rapidly fluctuating phase factor. These 
assumptions are similar to those made in the Chalker- 
Coddington model. For simplicity we suppose E = 0. ft 
is shown in Appendix A that this entails no loss of gen- 
erality or modification of behaviour; but the averaging is 
simpler to describe. 

To calculate the average Landauer conductance, we 
must average the S- matrix Ssusy(xf,Xi) over the ensem- 
ble described above. This is accomplished by expanding 
the S-matrix in a formal series (eq 28) and averaging term 
by term. The result is 

[Ssusvixf, Xi)] cns = exp[-(D/2)H C c{xf - Xi)] (48) 

where 



H cc = AB + BA. 



(49) 



A and B are defined below cq (42). The ensemble aver- 
aged conductance is then given by 



Mens = y < 0\c A c R exp f-^Hcc(xf - x t ) ) cf cf |() 



(50) 

Thus in order to calculate the average conductance we 
need to study the evolution of a two fermion state under 
the effective Hamiltonian Hcc- In contrast to Hsusy, 
Hcc is not random; but the price paid is that it is inter- 
acting. 



B. Anderson Localization 

The evaluation of eq (50) is simplified by the obser- 
vation that all of the fermionic bilinears that appear in 
Hcc (namely cf c Rj( , c R c R , cf c At and c A c A ) annhilate 
the two fermion state. Hence the expression for the av- 
erage conductance, eq (50), simplifies to 



[g]en S = < 0| B exp 



D 



rrboson,' 

H CC [ X f - X i 



|0 > B ■ 

(51) 



Here H c °^° n is the purely bosonic part of Hcc (written 
explicitly below) and |0 >b is the boson vacuum. Physi- 
cally then we need only calculate the boson vacuum am- 
plitude. 

Now let us analyse the boson Hamiltonian 



H 



cc 



= (&f b R] - b A b A )(b^ - b^b^ 1 ) 

-(b R b R - &f fe At )(6f - b A b^). (52) 



In contrast to the full Hamiltonian Hcc, the bosonic part 
is Hermitean. Note that it is also positive definite since it 
is of the form D + £>£>t wnere D = (b R b R - b^b Af ). 
After some manipulation we obtain the more revealing 
form 



rrboson 

H cc 



h n + h + + h-\ with 
h n = 2n R n R + 2n A n A + n R + n R + n A + n A + 2; 
h + = -2&f 6«t 6 A t ^t. 

h- = -2b R b R b A b A . 



(53) 



Here = fcf etc. Although lengthy eq (53) has 
a simple content. Consider the bosonic state \n >b, 
n = 0,1,2,... \n >b is defined as a normalized state 
that contains n bosons of each kind (R+, R-, A+ and 
A-). n — corresponds to the boson vacuum. Inspection 
of eq (53) shows that these states are closed under the 
action of H c °^ on . Hence we need only consider the block 
of -ffpc° n that lies within the invariant subspace spanned 
by these states. 

Our plan therefore is to find the eigenstates of H c °^° n 
that lie within the subspace spanned by \n >b- Ex- 
pansion of the boson vacuum in terms of these eigen- 
states will then allow straightforward evaluation of the 
boson vacuum amplitude and the average Landauer con- 
ductance, eq (51). Note that the Hermiticity of H c °S° n 
ensures that it's eigenstates form a complete set and it is 
therefore appropriate to use them as a basis. 

The coupled boson Hamiltonian, H c °^ on , is solved in 
section C below. It is found to have a continuum of eigen- 
states denoted \k > with eigenvalue (1 + k 2 )/2. Here 
k e [0, oo]. The eigenstates are orthogonal and arc nor- 
malized so that < k\0 >b= 1 and 



2 2 / irk 
< p\k >= — cosh 

7T 



2 J ksmh(^ 



-5(k-p). (54) 



It follows from eq (54) and the presumed completeness 
of the cigenfunctions of H c °, 



boson 
L CC 



tt [°° , fcsinhf 

- / dk n-A-\k >< k\ = X. 

2 J cosh 2 f 



(55) 



Here X denotes the identity matrix in the subspace 
spanned by \n >b- 

Inserting the resolution of the identity eq (55) into the 
conductance formula, eq (51), yields 



[fflons = exp y~^(xf ~ Xi) 

[°° k sinh ^ 
x / dk- 
Jo cosh 



T^k ex P 



Dk' 



-(x f - xA^j . (56) 



In the limit of large sample size we find 
[9Us " T { D(xf- Xi ) ) 2 ° XP {-J {Xf Xi) ) ■ (5?) 
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Thus the model exhibits Anderson localization: the con- 
ductance decays exponentially for sufficiently large sam- 
ple size as generally expected of a dirty one dimensional 
quantum wire. Eq (57) agrees with the result obtained 
byref §. 



C. Solution of Coupled Boson Hamiltonian 

We wish to solve the Schrodinger equation 

H%g*\4> >= A|V > (58) 

within the subspace spanned by \n >b- Expand the 
eigenstate as 

\tp >= a Q \0 > B +ai\l > B +■■ . + a n \n > B +■■■ (59) 

The effect of H^ on , eq (53), on the states \n >b is easily 
computed: 

h n \n> B = {An 2 +An + 2)\n> B ; 
h+\n> B = -2(n+ l) 2 \n + 1 > B ; 
h+\n > B = -2n 2 \n- 1 > B ■ (60) 

Eq (59) and (60) together allow us to write the 
Schrodinger equation as a three term recurrence relation 

(An 2 + An + 2)a n — 2n 2 a n -\ - 2(n + l) 2 a n +i = Aa„ 

(61) 

subject to a_i = 0. Our goal now is to solve eq (61) for 
different A and then ortho-normalize the solutions. The 
process of ortho-normalization will weed out the disal- 
lowed values of A. Note that from the Hermiticity and 
positive definiteness of Hqq 011 mentioned following eq 
(52), we are already assured that the allowed A must 
be positive and real. 

To solve the recurrence relation we introduce the gen- 
erating function 



f(x) = a + a\X + . . . + a n x" + 



(62) 



From the recurrence relation, eq (61), we can easily con- 
struct the differential equation obeyed by /: 

x ( x ~ !) 2 TT/ + ( Sx ~ !)(* - !)t-/ + ( x ~ 1 + f = 



This equation has three regular singular points: at x = 
0, 1 and oo; it is therefore a Riemann P-equation. One 
solution is analytic at x — as can be verified by directly 
substituting the series, eq (62), in the differential eq (63). 
This is the solution we seek; it generates the solution to 
eq(61). 

The solution to a Riemann P-cquation can always be 
expressed in terms of Hypergeometric functions. Making 
the standard transformations (see, e.g., ref |17|, chapter 
5) we find the analytic solution is given by 



f(x) = (l-x)"F(/i+l,/i + l,l;z); 



where jj, = — 



1 v 7 ! - 2A 



(64) 



2 2 

Note that for A > 1/2, /i becomes complex. Eq (64) 
therefore points to some change in behaviour at A = 1/2. 
Below we will find that in fact the allowed eigenvalues 
are A > 1/2. 

The coefficients a n can be extracted from f(x) via the 
contour integral 

- (65) 

Here C is a contour that encircles the origin but not the 
branch point x = 1. This contour integral cannot be 
expressed in elementary form in general but the large 
n behaviour of a n can be obtained by using an integral 
representation for the hypergeometric function (see Ap- 
pendix B). 

Let us focus on the solutions to eq (61) with A > 1/2. 
Write A = 1/2 + k 2 /2 where k <E [0, oo]. The large n 
asymptotic behaviour is 

1 

: COS I 

2 



2 , nk 
a n « - cosh — -pr —= 



(2fcsinh^) 



I ( k 

cos I — Inn + (j)(k) 



(66) 

Here the phase 4>(k) = arg[r(ifc)r(l/2 - ifc/2)/r 2 (l/2 + 
ik/2)]. Thus the solutions for A > 1/2 decay slowly ( as 
1/s/n ) and exhibit weak logarithmic oscillations. 

Fortunately, it turns out that this asymptotic be- 
haviour is sufficient to orthonormalize the eigenstates. 
Let b n be a solution to eq (61) with eigenvalue p: 

(An 2 + An + 2)b n - 2n 2 b n - X - 2(n + l) 2 b n+1 = pb n . 



dx 



dx' 



(63) 



(67) 

0. Multiply eq (61) by b n , eq (67) by a n , take the differ- 
ence and sum on n. A series of cancellations allows us to 
express the sum as a surface term 4 



4 This result is not surprising when one considers the analo- 
gous result of Strum-Liouville theory where it is known that 
the overlap integral of distinct eigenfunctions can be expressed 
as a surface term. See, for example, ref Il7|, p 719-720. 
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JV 



(P-A)£- 



2(N + l) 2 (a N+1 b N -a N b N+1 ). (68) 



n=0 



Thus the asymptotic behaviour of a n is sufficient to eval- 



uate the orthonormalization sum lim 
Use of eq (66), 
tion 5 



N- 



N 

n=0 u n u n . 



(68) and the delta function representa- 



lim — smkL = ird(k) 

L—>oo k 



(69) 



yields eq (54). 

We have so far focussed on solutions with A > 1/2. Us- 
ing similar arguments one can show that solutions with 
< A < 1/2 do not decay fast enough for large n to 
be ortho-normalizable even as continuum eigenfunctions. 
They are analogous to the exponential negative energy 
solutions of the free particle Schrodinger equation in el- 
ementary quantum mechanics. In the same spirit, the 
solutions for A > 1/2 are analogous to the positive en- 
ergy continuum of plane wave solutions. 

Having shown that there are no eigenstates with A < 
1/2, we have justified the step from eq (54) to eq (55). 
Alternatively one might try to prove the completeness re- 
lation, eq (55), directly. Forming the matrix element of 
eq (55) between the states \n >b and \m >b, one should 
check whether 



k sinh 

dk 



cosh' 



2 7rk 



a m (k)a n (k) = S„ 



(70) 



We have verified this analytically (and more extensively 
numerically) for a few small values of n and to for which 
the a's can be computed directly from eq (61). It may be 
possible to construct the general proof using eq (64-65), 
but this passes beyond conventional standards of rigour 
and good taste in theoretical physics. 



IV. DELOCALIZATION TRANSITION IN DYSON 
GLASS 



Now let us analyse the case in which m(x) is purely 
imaginary. We rewrite eq (2) as 



d ( $ 



' dx \ip- 



E —im \ / ?/M 
— im —E J \ ip- 



(71) 



Here m{x) is a real white noise process with mean mo. 
We have recycled the symbol m to conserve the finite re- 
sources of the alphabet. As noted in the introduction, 
this model is related to a model of glass first analysed by 
Dyson |L4| and to a special version of the random bond 
Ising Model introduced by McCoy and Wu |ll|. It has 
been extensively studied previously in its various incar- 
nations p^Jl^ l and it is known to have a critical point at 
rno = and E = 0. Here we shall use the Landauer con- 
ductance to investigate this critical point. For simplicity, 
the discussion is limited to two circumstances. First we 
set E — and study the behaviour of the conductance as 
too — > 0. Next we set too = and tune the Fermi energy, 
E. 



A. Elementary Solution for E = 

The techniques developed in this paper can be used to 
calculate the disorder averaged Landauer conductance of 
this model; but for the special case E = a much more 
complete solution can be obtained by elementary means. 
The content of this solution is instructive and it is there- 
fore described. 

When E = the transfer matrix of eq (71) has a par- 
ticularly simple explicit form for arbitrary m{x). Define 

M = I dxm(x). (72) 



Then the transfer matrix 

_ ( cosh M sinh M 
= \ sinhM cosh A/ 

Use of eq (4), (5) and (73) immediately yields 

t = sech M 

and 



.9 



—sech 2 
h 



M. 



(73) 



(74) 



(75) 



Since to is a Gaussian white noise process M is a Gaus- 
sian random variable. If we suppose 

Mens = m and 
[5m(x)5m(x')) cns = D6(x — x), (76) 

where 5m(x) = m(x) — mo, it is easy to calculate 

[M] ens = m (x f - Xi) and 



[{5Mf 



Dixi 



Xi) 



(77) 



From eq (77) the full distribution of M can be recon- 
structed: 



Used for example to obtain the orthonormalization of plane 
waves. 
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P(M) 



— Xi 



■. exp 



[M - mpjxf -Xj)X' 



2D(xf — Xi) 



(78) 



It is now altogether straightforward to obtain statisti- 
cal information about the Landauer conductance from eq 
(75) and (78). In particular it can be shown: (i) At the 
critical point the average conductance decays as a power 
of the system size (xt — x^ 1 ! 2 . (ii) Away from the 
critical point (m =/= 0) the average conductance decays 
exponentially (Anderson localization). The localization 
length diverges as m^ v where the exponent v = 2. (hi) 
The conductance has a very broad distribution; the av- 
erage is dominated by large rare fluctuations, (iv) Away 
from the critical point the typical conductance also de- 
cays exponentially but with a distinct localization length 
that diverges with exponent v = 1 as the critical point is 
approached. 

In summary although eq (71) is essentially trivial to 
solve for E = 0, the content of the solution is not triv- 
ial. It illustrates many general features of random critical 
systems Q. 



B. Tuning Fermi Energy 

1. Disorder Averaging 



\t\ 2 =< 0\c*S§(xf,Xi)c™\0 X 0\b^S§(x f ,Xi)b^\0 > . 

(79) 

Having noted these possibilities, we now disregard 
them. Instead we use eq (44) to calculate \t\ 2 . Hence 
we need to average the S-matrix Ssusv(xf, Xi) over the 
white noise ensemble of real m(x) with [m(a;)] cns = and 
[m(x),m(x')] ens = DS(x — x'). The result is 



[SsusY{xf,Xi)] ens = exp -—H D (xf - x t ) 



D 



where 



jN- 



D 



K 



(80) 



(81) 



N counts the total number of particles regardless of 



- At «t «t 



C R C R + C A\ C A^ + C A C A _ 



species and K = i(c, 
6+ f &5 + b R bf - 6+ V f + b A bf). The ensemble averaged 
conductance is then given by 



[g] ens = < 0|c+c+exp 



D 

—H D (x f - Xi) 



-+ c + 



|Q> 
(82) 



The localization length is expected to diverge as 
(In E) 2 as E — * 0. To probe this singularity we shall make 
one simplification. The frequency E will be considered to 
be complex, and we shall approach the singularity at the 
origin along the imaginary axis. Thus we shall replace 
E — * iuj in eq (71) and study [ff] e ns(w). 

Note that for imaginary frequency eq (71) is real and 
therefore the scattering amplitude t is also real. In gen- 
eral t is not a suitable probe of delocalized behaviour as 
it decays exponentially (due to phase fluctuations) even 
when \t\ 2 does not. In this case however we may look 
for the transition by calculating the average of t which is 
not only real, but also positive (and therefore free of sign 
fluctuations 6 ). This constitutes a simplification because 
to calculate t the A fermions and bosons are not needed. 

Even if one insists on calculating \t\ 2 , the A fermions 
and bosons are unneeded. Since t is real \t\ 2 = t 2 and we 
may write 



2. Critical Behaviour 

Once again all the fermionic bilinears in K 2 annhilate 
c+ c+ |0 >; hence the conductance is determined by the 
boson vacuum amplitude 

e 2 

[g]ens = -T exp[-2w(x/ - Xi)] 

x < 0| fl exp (-^H h D os ™{x f - Xi )j |0 > B . (83) 

As before H]^ son is hermitean and positive definite. It 
has an invariant subspace spanned by the states |n >b, 

n = 0, 1, 2, \n >b is here defined as the state with 2 n 

bosons of each kind. Hence we may focus on the block of 
-ffboson that lies within the subspace spanned by \n >b- 



Proof: Consider a series of truncated problems as in section 
II C. First note that for xo = x/, t(xo) = 1. Then observe 
that t(xo) can never vanish; if it did, the wave function van- 
ishes in the drain, and by integrating eq (71) backwards, it 
vanishes in the sample and source also. Since it cannot pass 
through zero, by continuity, t(xo) must remain positive for all 
x including x t . 
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Our plan is to find eigenstates of this block and to ex- 
pand the boson vacuum in terms of these eigenstates. In 
contrast to the two edge Chalker-Coddington case the 
eigenstates are discrete and can be labelled by an in- 
teger I = 1,2,3,.... The eigenvalue problem is solved 
approximately in the relevant small u> limit in the next 
subsubsection. It will be found that the eigenvalues are 
given by 



2/2 



A, 



2tt 2 Z 

WW 



and the overlap with the vacuum is 

1 sinh nki 



<i\o> B r = 



2 In 



D 



nk. 



(84) 



(85) 



where ki = lir / \a(u) / D) . Putting together these results 
e 2 1 ( _ o„ 1 



h 12 In ■ 



exp \-2tt z D 



Xi) 



(86) 



for large samples. Eq (86) shows that the sample is An- 
derson localized away from the critical point. The lo- 
calization length diverges as (lna>) 2 as u> — > 0, in agree- 
ment with previous work (Balents and Fisher |fl3f and 
refs therein). 



3. Solution of Boson Problem 

This calculation is very similar to the corresponding 
calculation by Balents and Fisher |T^] of the average 
Green's function . For this reason, after pointing out the 
special feature of the conductance calculation (discussion 
leading to eq (89)) we present the remaining steps in out- 
line. The reader interested in more details should consult 
the appropriate sections of Balents and Fisher. 

We wish to solve 



However if we work with an unusual normalization 
such that < n|s|m >b~ (n+ 1)6 mn simple coefficients 
result. Adopting this convention we find 



~~D 



m + 4m + 4m + 1 a 



— m(2m — l)a m _i — (m + l)(2m + 3)a m+ i = Xa m . (89) 

The corresponding generating function obeys a differ- 
ential equation with three singular points, two regular 
and one irregular. Thus the solutions are related not 
to the common functions of mathematical physics but 
to the more obscure Mathieu or Spheroidal functions. 
Rather than pursue this direction we closely follow Ba- 
lents and Fisher to obtain an approximate solution to the 
Schrodinger eq (89), valid in the interesting limit of small 
lo/D. 

First we solve eq (89) for to = using the generating 
function method. This solution should be accurate for 
n <C D/ui. Setting clq — 1 we find 



1 / sinh irk 
2n \ irk 



1/2 



cos[fc Inn — 4>(k)} 



(90) 



for 1 <C n <C D/uj. Here we have introduced k via 
A = 2fc 2 and 4>(k) = argT(l - ik) + 2fcln2. Next we 
approximate the large n behaviour by taking the contin- 
uum limit of eq (89): 



* d 
on— a - 

an 



D 



(91) 



Eq (91) is soluble by Laplace transforms; but an even 
quicker solution may be effected by introducing f(n) = 
na(n) which obeys the equation 



rrboson 



\ip >= X\ip > 



(87) 



within the subspace spanned by the states \n >b defined 
in the previous subsubsection. To this end expand \ip > 
as 



\ip >= a \0 > B +<2i|l > B + ■ ■ ■ + a n \n > B + . 



(88) 



A difficulty arises if we assume the states are normal- 
ized so that < n|s|m >b— S mn . In this case the effect 
of i/^, oson on | n >b is to yield a linear combination of 
\n + 1 >Bi \ n >b aud \n — 1 >b with coefficients that 
involve products such as \Jn{n + 1). Such non-integer 
coefficients would then also appear when eq (87) is writ- 
ten as a recurrence relation analogous to eq (61) and 
would defeat the generating function technique that was 
used to solve eq (61). 



n 2 — /• 
dn 2 



n—f+(---n)f = 0. 



dn" 



D 



(92) 



Precisely this equation was analysed by Balents and 
Fisher in Appendix D of their paper. Borrowing their 
results we can find the values of A for which eq (91) has 
a solution that decays as n — > oo and smoothly matches 
eq (90) for 1 <C n <C D/lj. This yields the quantized 
energy levels of eq (84). 

The eigenstates may be approximately normalized by 
taking the continuum solution to apply everywhere. The 
normalization sum X)^°=o a n/( n +l) ma y then be replaced 
by an integral. The relevant integral is evaluated asymp- 
totically by Balents and Fisher (Appendix D). Transcrip- 
tion of their result yields eq (85). 
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C. Summary 



In this interaction representation it is easy to show 



For reference we summarize the results of this section: 
At the critical point (mo = 0, to = 0) the average con- 
ductance decays as the inverse square root of the sample 
size. Off the critical point the average conductance de- 
cays exponentially on a length scale called the localiza- 
tion length. For w = wc find the localization length 
of the average conductance diverges as rn^ u as mo — > 0, 
with exponent v = 2; for mo = 0, we find it diverges as 
(lnw) 2 . These results for the average conductance local- 
ization length are consistent with previous calculations 
of other non-local correlation functions for these models 

IIS 

For u> = 0, the model is completely soluble by ele- 
mentary methods. In this simple case it shows many of 
the features generally expected of random critical sys- 
tems In particular the Landauer conductance is 
very broadly distributed and the average is dominated by 
large, rare fluctuations. Away from the critical point the 
typical conductance also decays exponentially but with 
a localization length different from the average conduc- 
tance. The localization length for the typical conduc- 
tance diverges with exponent 1 as too — > 0. 
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A[x) = exp[-i2E(x - Xi)]A, 
B(x) = exp[i2E{x - Xi)]B. 



(A4) 



The disorder average of U can be computed as in sec- 
tion III. 

[U{x,Xi)] eas = exp ( - — (x - x t )(AB + BA) J . (A5) 



Recall AB + BA = H cc (eq 49). Hence 



[SsusY{x,Xi)] ens = exp[iEM(x - X*)] 

D 



xexp[-—Hcc{x-x l )). (A6) 



Finally note 

Mcfc^\0>=0. (A7) 
Eq (44), (A6) and (A7) together show that 



Mens = \< 0|c+4'exp (-£-H CC (Xf - Xi) ) cfcf\() ■ 



(A8) 



-independent of E as claimed in the paper. 
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APPENDIX A: 

The purpose of this Appendix is to show that the re- 
sults of section III are not modified when E ^ 0. This 
can be done in various ways. The method followed here 
is chosen because it also illustrates how to average the 
S-matrix when the random Hamiltonian, H$usy, is com- 
posed of two non-commuting pieces. 

Factorise the S-matrix as 

SsusY{x,Xi) = exp[iEM(x - Xi)]U(x,Xi). (Al) 
Then eq (42-43) imply that U obeys 
d 

- i—U{x,Xi) = [m(x)A(x) + m*(x)B]U(x,Xi) (A2) 

and the initial condition U{xi,Xi) — 1. Here 

A(x) = exg>[-iEM{x - Xi)]Aexp[iEM(x - a;*)], 

B(x) = exp[-iEM(x - Xl )]B exp[iEM(x - a*)]- (A3) 



APPENDIX B: 

The contour integral in eq (65) is evaluated here 
asymptotically for large n. /(x) is defined in eq (64). 
As in the paper, we focus on A > 1/2 and write A = 
(1 + fc 2 )/2. Thus fj, = -1/2 + ik/2. 

Use the integral representation of the hypergeometric 
function 



F(a, b, c; x) 



r(c) 



r(&)r( c -6) A 



dt{t-x)~ a t a - c {t- 1) 



c-b-1 



(Bl) 



valid provided Rec > Reb > (ref [Q, chapter 5). 

Substitute eq (Bl) in eq (65) with a — > /U+l, b — > fi+X 
and c — ► 1. Exchange the order of the t and x integrals 
to obtain 



r(c) 



r(&)r( c -&) 



dte-^t-if-^ut) (B2) 



where 



/»(*) 



dx 1 



2iri x n+ 



(B3) 



Here C is a contour that encircles the origin but not the 
branch point at x = 1. 
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Note that the integrand in eq (B3) has branch points 
at x = 1 and x = t where t is some point on the positive 
real axis to the right of 1 (and to be eventually integrated 
over the range from 1 to oo). To be consistent with the 
conventions of the integral representation, eq (Bl), we 
draw branch cuts along the positive real axis from 1 to 
+oo and t to +oo. Also the phase of (t — x) and (1 — x) 
must both be taken to be zero when x lies on the real 
axis to the left of 1. 

The contour C is now deformed to pass above and be- 
low the branch cut. It is closed by a small circle around 
x = 1 and by a big circle at oo. The contribution of the 
circles to the contour integral vanishes. The contribution 
from integrating above and below the branch cut is 



F{a, b, c; x) 



T(c)T(c- 



b) 



T(c-a)T{c-b) 
T(c)T{a + b- c) 



F(a,b, a + b- c + l;l - x) 
(l-x) c -°- b 



T(a)T(b) 

xF(c- a,c- b,c- a - b+ 1; 1 - x). (B8) 



We obtain 

,1 k 1 
F(-+i-, - 
V 2 2' 2 



k 

• i— , 1: 
2' ' 



r(-ijfc) 



r 2 (i/2-ifc/2) 
r(ifc) 

+ THl/2Tik/^j ! 



(B9) 



„ . . 1 , 7rfc 
f n (t) = - cosh — 

7T Z 



dx- 



1 



„n+ 



-l/2-ifc/2 



(x-1) 



-l/2+ifc/2 



Substitute eq (B9) in eq (B7). 
eq (66). 



This leads to the result of 



(B4) 

Once again exchange the order of the x and t integrals 
to obtain 



1 



cosh' 



irk 



i 



dtt-Va+<fc/2( t _ _ x yi/2- ik /2_ (B5) 



The r function prefactors in the integral representation, 
eq (Bl), have been simplified by use of the formula 
r(z)r(l -z) = 7r/sin(7rz). 

Rescale the t integral so the lower limit becomes 1. 
Upon comparison with the integral representation, eq 
(Bl), it is seen 

1 

dx — —(x - iy 
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